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Abstract. We present the redshift distribution of a complete sample of 480 galaxies with K s < 20 distributed over 
two independent fields covering a total area of 52 arcmin . The redshift completeness is 87% and 98% respectively 
with spectroscopic and high-quality and tested photometric redshifts. The redshift distribution of field galaxies 
has a median redshift 2 me d ~ 0.80, with ~ 32% and ~9% of galaxies at z > 1 and z > 1.5 respectively. A 
"blind" comparison is made with the predictions of a set of the most recent ACDM hierarchical merging and pure 
luminosity evolution (PLE) models. The hierarchical merging models overpredict and underpredict the number 
of galaxies at low-z and high-z respectively, whereas the PLE models match the median redshift and the low- 2 
distribution, still being able to follow the high-z tail of N(z). We briefly discuss the implications of this comparison 
and the possible origins of the observed discrepancies. We make the redshift distribution publicly available. 
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1. Introduction 

The mass assembly history of galaxies remains one of 
the critical issues in observational cosmology: did galax- 
ies reach their present stellar mass only recently (say, at 
Z ,$ 1) ? Or were most (massive) galaxies already in place 
by z ~ 1 ? Spectroscopic surveys of faint galaxies selected 
in the _ftT-band currently offer the best opportunity to an- 
swer these questions (Broadhurst et al. 1992). The main 
advantages with respect to optically selected samples in- 
clude: the direct sensitivity to the galaxy stellar mass 
rather than to the ongoing/recent star formation activ- 
ity (Gavazzi et al. 1996; Madau et al. 1998), the smaller 
K-correction effects, and the minor influence of dust ex- 
tinction. 

In this framework, we have completed an optical 
and near-infrared spectroscopic survey down to K s < 
20 (dubbed "K20 survey") using ESO VLT telescopes 
and instruments, with full survey details being given 
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in Cimatti et al. (2002b, hereafter Paper III; see also 
tittp : //www. arcetri . astro . it/~k20/). The K20 sam- 



ple includes 546 objects to K s < 20 (Vega system) over 
two independent fields (52 arcmin 2 in total), so to be 
less affected by the cosmic variance. The spectroscopic 
redshift completeness is 94% and 87% for K s < 19 and 
K s < 20 respectively. This makes the K20 sample the 
largest and most complete spectroscopic sample of galax- 
ies with K s < 20 available to date (see Paper III; cf. Cowie 
et al. 1996; Cohen et al. 1999). Moreover, a 98% redshift 
completeness is reached for the K s < 20 sample when in- 
cluding the photometric redshifts obtained with the avail- 
able deep UBVRIzJKg imaging for those objects without 
a spectroscopic redshift. If stars and broad- line AGNs are 
excluded, the total number of galaxies with K s < 20.0 and 
with redshifts is 480. 

In two previous papers based on the K20 survey we 
showed that Extremely Red Objects (EROs, defined by 
R — K s > 5) are nearly equally populated by old passively 
evolving galaxies and by dusty star-forming systems at 
z ~ 1 (Cimatti et al. 2002a, Paper I). The number of all 
(old+dusty) EROs is strongly underpredicted by hierar- 
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chical merging models (HMMs), whereas old EROs have a 
density consistent with PLE models for passive early-type 
galaxies (Paper I), and are strongly clustered as opposed 
to dusty EROs (Daddi et al. 2002, Paper II). 

In this Letter, we present and discuss the observed 
redshift distribution, N(z), of all the galaxies in the K20 
sample, irrespective of their color, and compare it to the 
expectations for the case of PLE of galaxies, as well as to 
the predictions of various HMM renditions. The currently 
favored cosmological model is adopted, i.e., Hq = 70 km 
s- 1 Mpc" 1 , n m = 0.3 and fl A = 0.7. 

2. Observations vs. model predictions 

The observed N(z) for the 480 galaxies (417 with 
spectroscopic and 63 with photometric redshifts 
respectively) with K s < 20 is shown in Fig. 1. 
The redshift distribution can be retrieved from 
http : //www, arcetri . astro . it^k20/releases| The 
spike at z ~ 0.7 is due to two clusters (or rich groups) 
of galaxies respectively at 0.665 < z < 0.672 (23 galax- 
ies) and 0.732 < z < 0.740 (33 galaxies) (see Paper 
III). The median redshift of N(z) is z me d — 0.737 
and z me d — 0.805, respectively with and without the 
two clusters being included. Without the clusters, the 
fractions of galaxies at z > 1 and z > 1.5 are 138/424 
(32.5%) and 39/424 (9.2%) respectively. The high-z tail 
extends beyond z = 2. The contribution of objects with 
only a photometric redshift becomes relevant only for 
z > 1.5. The fractional cumulative distributions displayed 
in Fig. 2-3 (bottom panels) were obtained by removing 
the two clusters mentioned above in order to perform 
a meaningful comparison with the galaxy formation 
models which do not include clusters (PLE models) or 
are averaged over very large volumes, hence diluting the 
effects of redshift spikes (HMMs). 

No best tuning of the models was attempted in this 
comparison, thus allowing an unbiased "blind" test with 
the K20 observational data. The model predicted N(z) 
have been normalized to the K20 survey sky area. We 
first discuss the case of PLE expectations, as derived by 
Pozzctti et al. (1996, 1998, PPLE hereafter), and Totani 
et al. (2001) (TPLE hereafter). 

2.1. Comparison with PLE models 

In the PPLE model the present-day galaxy luminosity 
function is divided into 5 Hubble types (E, SO, Sab- 
Sbc, Scd-Sdm, Im). The spectral evolution for each type 
is described by the Bruzual & Chariot (1993) model 
(GISSEL version 2000) which reproduces the rest-frame 
colors and K-corrections of local galaxies. Exponentially 
declining star formation rate histories and solar metallic- 
ity are adopted. The age of each galaxy is set to 12.5 Gyr 
(zf = 5.7) with the exception of Im galaxies (age=0.1 
Gyr). The e-folding times are set to 0.3, 2, 10 Gyr, and 
oo for E, SO, Sab-Sbc, and Scd-Sdm-Im galaxies respec- 
tively. Dust extinction is not taken into account. Besides 
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Fig. 1. The observed differential redshift distribution. The 
shaded histogram shows the contribution of photometric red- 
shifts. The bin at z < indicates the 9 objects without redshift. 

the adopted cosmology, the only difference between the 
PPLE model used here and Pozzetti et al. (1996, 1998) 
is the use of the if s -band local luminosity function from 
the 2MASS for different morphological types (Kochanek 
et al. 2001), which is in agreement with the overall local 
luminosity function of Cole et al. (2001). Fig. 2 shows the 
predictions of the PPLE models for two types of initial 
mass function (IMF), Salpeter (1955) and Scalo (1986). 
With the flatter (i.e. Salpeter) IMF the intrinsic luminos- 
ity of both passively evolving and star forming galaxies 
increases more rapidly with redshift than in the case of 
the steeper (Scalo) IMF. As already discussed by Pozzetti 
et al. (1996) (see also McCracken et al. 2000), the mild evo- 
lution allowed by the Scalo IMF is more consistent with 
the observations of the rest-frame ultraviolet luminosity 
density up to z ~ 1 (e.g. Cowie et al. 1999), and repro- 
duces most observables (galaxy counts, color and redshift 
distribution from the U to the K band) without invoking 
the strong number density evolution or dust extinction 
required by the Salpeter IMF (cf. Fig. 2). 

In the TPLE model galaxies are also divided into 5 
types (E/S0, Sab, Sbc, Scd, Sdm) and their spectral evolu- 
tion is described using the Arimoto & Yoshii (1987) mod- 
els. The £?-band local luminosity function is used, and the 
Salpeter IMF and a top-heavy IMF with exponent 0.95 are 
adopted for spirals and ellipticals respectively, with two 
options for formation redshifts: zj=3 and z/=5. Contrary 
to the PPLE model, the evolution of metallicity and dust 
extinction in galaxies is treated in this model. 

Fig. 2 shows fairly good agreement between the ob- 
served N(z) distribution and the PLE models (with the 
exception of PPLE with Salpeter IMF). The predicted 
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Fig. 2. Top panels: the observed differential N(z) for K s < 
20 (histogram) compared with the PLE model predictions. 
Bottom panels: the observed fractional cumulative redshift dis- 
tribution (continuous line) compared with the same models. 
The left and right panels show the models without and with 
the inclusion of the photometric selection effects respectively. 
Sc and Sp indicate Scalo and Salpeter IMFs respectively. 



and the observed total number of galaxies with K s < 20 
agree within 10% for the PPLE and the TPLE (z f = 3), 
28% for the TPLE (z f = 5), and 34% for PPLE with 
a Salpeter IMF. The predicted median redshifts are just 
slightly higher than observed: z me d— 0.83, 0.75 and 0.79 for 
the PPLE, TPLE (zf = 3, 5) models respectively, but in- 
consistent with the PPLE with Salpeter IMF (z med —1.05). 
This is due to these PLE models somewhat overpredicting 
the number of galaxies at z ^ 1.2. However, as extensively 
discussed in Paper III, because of the photometric selec- 
tion effects present in the K20 sample (partly due to the 
cosmological surface brightness dimming), the total fluxes 
of spirals and ellipticals with L ~ L* (i.e. the bulk of the 
K20 sample) are, on average, underestimated by about 0.1 
and 0.25 magnitudes, respectively. In order to assess the 
influence of such effects, we compared the observed red- 
shift distribution (down to our nominal K s < 20.0) with 
the PPLE and TPLE models with K s < 19.9 for "disk" 
and K s < 19.75 for "early-type" galaxies. Fig. 2 (right 
panels) shows that when such selection effects are taken 
into account the PLE models become much closer to the 
observed N(z) thanks to the decrease of the predicted 
high-z tail. According to the Kolmogorov-Smirnov test, 
the PLE models are acceptable at 95% confidence level, 
with the exception of the PPLE model with Salpeter IMF 
(rejected at > 99% level). We conclude that PLE models 
offer a satisfactory fit to the observed N(z) distribution, 
all the way to the highest redshifts in our sample. 



2.2. Comparison with hierarchical merging models 

For the comparison with the HMM predictions we were 
kindly granted access to the model databases of Cole at 
al. (2000, COO hereafter), Somerville et al. (2001, S01 here- 
after) and Mcnci et al. (2002, M02 hereafter). These mod- 
els are tuned to reproduce some low- z observable, such as 
the local galaxy luminosity function near L ~ L* (COO, 
M02) or the Tully-Fishcr relation (S01). The main differ- 
ence among the HMMs used here is the inclusion in S01 
of the merging-promoted "starburst" mode of star forma- 
tion besides the "quiescent" mode, the only one included 
in COO and M02. The starburst mode has the effect of in- 
creasing the overall star formation at high redshift, when 
most of the merging takes place. Moreover, merging be- 
tween satellite galaxies within DM haloes is included in 
S01 and M02, but neglected in the COO, where satellites 
are allowed to merge only onto the central massive galaxy. 
The M02 model without merging between satellites is oth- 
erwise equivalent to the COO rendition. The effect of merg- 
ing between satellites is to deplete the number of low-mass 
galaxies which aggregate to form larger units, thus flatten- 
ing the galaxy mass function at the faint end and slightly 
increasing the number of intermediate mass galaxies (see 
S01 and M02 for more details). All the HMMs used here 
adopt a Salpeter IMF. 

The HMMs overpredict the total number of galaxies 
with K s < 20 in the K20 survey area by factors of 30-45%. 
In particular, Fig. 3 (top panels) shows that all the HMMs 
show an excess of predicted galaxies at z < 0.5, e.g., by a 
factor of ~ 2.5 at z ~ 0.4 for the COO model, and - 1.5 - 2 
for the M02 and S01 models, respectively. The predicted 
median redshifts are z me d— 0.59, 0.70 and 0.67 for the COO, 
M02 and S01 models, respectively, thus being systemat- 
ically lower than the observed z me d- Moreover, all the 
HMMs have a deficit of z > 1 galaxies, in particular with 
a fraction at z > 1.5 smaller by factors of ~4 for the COO 
model and ^2-^3 for the M02 and S01 models. Such a dis- 
crepancy increases dramatically for higher redshifts, where 
all the HMMs predict no galaxies with z > 2. Fig. 3 (bot- 
tom panels) illustrates that in the fractional cumulative 
distributions the discrepancy with observations appears 
systematic at all redshifts. The Kolmogorov-Smirnov test 
shows that all the HMMs are discrepant with the obser- 
vations at > 99% level. The inclusion of the photometric 
biases exacerbates this discrepancy, as shown in Fig. 3 
(right panels) for the M02 model calculated for K s < 19.8 
in order to include an average photometric bias for spi- 
rals and ellipticals (the discrepancy for the COO and S01 
models becomes even stronger). 

The excess of galaxies at z 0.5 seen in Fig. 3 is due 
to HMMs predicting too many low-mass, low- luminosity 
galaxies. This excess has typically afflicted HMMs, with 
the merging between satellites improvement being appar- 
ently insufficient to provide a better agreement with the 
data. But in addition, HMMs underpredict the number of 
high-rcdshift objects. This is illustrated by Fig. 4, where 
the PPLE model is capable to reproduce the cumulative 
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Fig. 3. Top panels: the observed differential redshift distribu- 
tion for K s < 20 (histogram) compared with the HMM predic- 
tions. Bottom panels: the observed fractional cumulative red- 
shift distribution (continuous line) compared with the same 
models of top panels. The right panels show the M02 model 
with the inclusion of the photometric selection effects. 




Fig. 4. The observed cumulative number of galaxies between 
f < 2 < 3 (continuous line) and the corresponding poissonian 
±3cr confidence region (dotted lines). The PPLE (Scalo IMF) 
and the M02 models are corrected for the photometric biases. 



number distribution of galaxies at 1 < z < 3 within l-2cr, 
whereas the M02 model is always discrepant at ~ 3<r level 
(up to > 5a for 1.5 < z < 2.5). All the results described in 
this section remain valid if shallower limiting magnitude 
thresholds are adopted (see Fig. 5). 

3. Discussion 

Early predictions of the expected fraction of galaxies at 
z > 1 in a K < 20 sample indicated respectively « 60% 
and « 10% for a PLE case and for a (then) standard 
Q m = 1 CDM model (Kauffmann & Chariot 1998). This 
version of PLE was then ruled out by Fontana et al. 
(1999). The more recent PLE and HMM used in this pa- 
per consistently show that for z > 1 the difference be- 
tween the predictions of different scenarios is much less 
extreme. These results partly from the now favored ACDM 
cosmology which pushes most of the merging activity in 
hierarchical models at earlier times compared to tCDM 
and SCDM models with fl m = 1, and partly to different 
recipes for the star formation modes, which tend to narrow 
the gap between HMMs and the PLE case (e.g. Somerville 
et al. 2001; Firth et al. 2002). The disagreement between 
the observed N(z) and the predictions of the most up- 
dated HMMs based on a ACDM cosmology would then 
become even stronger in the case of old-fashioned CDM 
models with fl m = 1 because structures form later in a 
matter-dominated universe, and thus they would predict 
an even lower fraction of galaxies at high-z. In this re- 
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Fig. 5. The observed differential redshift distribution for K s < 
19 and K s < 19.5 (histograms) compared with the model pre- 
dictions (not corrected for photometric selection effects). 



spect, our results can be seen as additional evidence that 
the universe is not matter-dominated (f2 m < 1), and sug- 
gest that the HMMs may perform better if Q m is even 
lower than the currently favored fi m = 0.3. 

Nevertheless, the results of the K20 survey indicate 
that the shape and the median of the observed redshift 
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distribution of K s < 20 galaxies are in broad agreement 
with the expectations of PLE models, while disagree with 
the predictions of current hierarchical merging models of 
galaxy formation. This discrepancy refers to all galaxies, 
irrespective of color or morphology selection, and therefore 
is more general than the already noted discrepancies with 
EROs (Daddi et al. 2000; Paper I; Cimatti 2002). The poor 
performance of HMMs in accounting for the properties of 
even z = — >~ 1 early-type galaxies has been emphasized 
in the past (e.g., Renzini 1999; Renzini & Cimatti 1999). 
Moreover, among low-redshift galaxies there appears to 
be a clear anti-correlation of the specific star formation 
rate with galactic mass (Gavazzi et al. 1996; Boselli et 
al. 2001), the most massive galaxies being "old", the low- 
mass galaxies being instead dominated by young stellar 
populations. This is just the opposite than expected in 
the traditional hierarchical merging scenario, where the 
most massive galaxies are the last to form. 

On the other hand, the strong clustering of EROs 
seems to be rather consistent with the predictions of CDM 
models of large scale structure evolution (Daddi et al. 

2001, Paper II; Firth et al. 2002). Thus, adopting the hier- 
archical merging ACDM scenario as the basic framework 
for structure and galaxy formation, the observed discrep- 
ancies may be ascribed to the heuristic algorithms adopted 
for the star formation processes and their feedback, both 
within individual galaxies and in their environment. Our 
results suggest that HMMs should have galaxy formation 
in a CDM dominated universe to closely mimic the old- 
fashioned monolithic collapse scenario. This requires to 
enhance merging and star formation in massive haloes 
at high redshift (say, z <L 3), while in the meantime sup- 
pressing star formation in low-mass haloes. For instance, 
Granato et al. (2001) suggested the strong UV radiation 
feedback from the AGN activity during the era of super- 
massive black hole formation to be responsible for the 
suppression of star formation in low-mass haloes, hence 
imprinting a "anti-hierarchical" behavior in the baryonic 
component. The same effect may well result from the feed- 
back by the starburst activity itself (see also Ferguson & 
Babul 1998). 

In summary, the redshift distribution presented in this 
paper, together with the space density, nature, and clus- 
tering properties of the ERO population (Paper I, Paper 
II) and the redshift evolution of the luminosity and stellar 
mass functions derived for the K20 sample (Pozzetti et al. 

2002, Fontana et al. 2002) provide a new set of observ- 
ables on the galaxy population in the z ~ 1 — 2 universe, 
thus bridging the properties of z ~ galaxies with those 
of Lyman-break and submm/mm-selected galaxies at z 
~ 2-3. While making a step towards the fully empirical 
mapping of galaxy formation and evolution, this set of 
observables poses a new challenge for theoretical models 
to properly reproduce. 

Acknowledgements. We are in debt with Carlton Baugh, 
Rachel Somerville and Tomonori Totani for providing their 
"blind" model predictions. We thank the referee, Nathan 



Roche, for useful comments and the VLT support astronomers 
for their assistance during the observations. AC warmly thanks 
ESO (Garching) for the hospitality during his visits. 

References 

Arimoto N., Yoshii Y. 1987, A&A, 173, 23 
Baugh CM., Cole S. & Frenk C.S., 1996, MNRAS 283, 1361 
Boselli A., Gavazzi G., Donas J. & Scodeggio M. 2001, AJ, 
121, 753 

Broadhurst, T., Ellis, R.S., & Grazebrook, K. 1992, Nature, 
355, 55 

Bruzual G., Chariot S. 1993, 405, 538 

Cimatti A. 2002, in "The Mass of Galaxies at Low and High 
Redshift", ESO/ USM Workshop Venice, Italy, Springer- 
Verlag, in press, istro-ph/020105C 

Cimatti A., Daddi E., Mignoli M. et al. 2002a, A&A, 381, L68 
(Paper I) 

Cim atti A., Mignoli M ., Daddi E. et al. 2002b, A&A, in press 

flastro-ph/0206168; ) (Paper III) 
Cohen J.G., Blandford R, Hogg DW. et al. 1999, ApJ, 512, 

30 

Cole S., Lacey C.G., Baugh CM. & Frenk C.S. 2000, MNRAS, 
319, 168 

Cole S., Norberg P., Baugh CM. et al. 2001, MNRAS, 326, 
255 

Cowie, L.L., Songaila A., Hu E.M. & Cohen J.G. 1996, 
AJ,112,839 

Cowie, L.L., Songaila A. & Barger A.J. 1999, AJ, 118, 603 
Daddi E., Cimatti A. & Renzini A., 2000, A&A 362, L45 
Daddi E., Broadhurst T., Zamorani G. et al. 2001, A&A, 
376,825 

Daddi E., Cimatti A., Broadhurst T. et al. 2002, A&A, 384, 

LI (Paper II) 
Ferguson H.C & Babul A. 1998, MNRAS, 296, 585 
Firth A.E., Somerville R.S., McMahon R.G. et al. 2002, 

MNRAS, 332, 617 
Fontana A., Menci N., D'Odorico S. et al. 1999, MNRAS, 310, 

L27 

Fontana A. et al. 2002, in preparation 
Gavazzi C, Pierini D. & Boselli A. 1996, A&A, 312, 397 
Granato G.L., Silva L., Monaco P. et al. 2001, MNRAS, 324, 
757 

Kauffmann C, Chariot S. 1998, MNRAS, 297, L23 
Kochanek C.S., Pahre M.A., Falco E.E. et al. 2001,ApJ,560,566 
Madau P., Pozzetti L. & Dickinson M., 1998, ApJ, 498, 106 
Menci N., Cavaliere A . , Fontana A . , Gial longo E. & Poli F. 

2002, ApJ, in press ( |a,stro-ph/0204178[ ) 
McCracken H.J., Metcalfe N., Shanks T. et al. 2000, MNRAS, 

311, 707 

Pozzetti L. Bruzual A.G. & Zamorani, G. 1996, MNRAS, 281, 
953 

Pozzetti L., Madau P., Zamorani G., Ferguson H.C. & Bruzual 

A.G. 1998, MNRAS, 298, 1133 
Pozzetti, L., et al. 2002 (in preparation) 

Renzini A. 1999, in The Formation of Galactic Bulges, ed. 
CM. Carollo, H.C. Ferguson, & R.F.G. Wyse (Cambridge: 
CUP), p. 9 

Renzini A., Cimatti A. 1999, ASP Conf. Ser. 193, 312 

Salpeter E.E. 1955, ApJ, 121, 161 

Scalo J.M. 1986, Fundam. Cosm. Phys., 11, 1 

Somerville R.S, Primack JR. & Faber S.M. 

2001,MNRAS,320,504 
Totani T., Yoshii Y., Maihara T., Iwamuro F. & Motohara K. 

2001, ApJ, 559, 592 



